Feshbach-Einstein condensates 
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We investigate the phase diagram of a two-species Bose-Hubbard model describing atoms and 
molecules on a lattice, interacting via a Feshbach resonance. We identify a region where the sys- 
tem exhibits an exotic super-Mott phase and regions with phases characterized by atomic and/or 
molecular condensates. Our approach is based on a recently developed exact quantum Monte Carlo 
algorithm: the Stochastic Green Function algorithm with tunable directionality. We confirm some 
of the results predicted by mean-field studies, but we also find disagreement with these studies. In 
particular, we find a phase with an atomic but no molecular condensate, which is missing in all 
mean-field phase diagrams. 
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More than eighty years ago, Einstein predicted a re- 
markable phenomenon to occur in a gas of identical 
atoms interacting weakly at low temperature and high 
density [ij. Under such conditons, when the de Broglie 
wavelength of the atoms becomes larger than the inter- 
atomic distance, a macroscopic fraction of the atoms ac- 
cumulates in the lowest energy state. This phenomenon, 
known as Bose-Einstein condensation, remained in the 
archives for a long time, and was reconsidered later with 
the discovery of the superfluidity of Helium in 1937. It 
is only in 1995 with the advent of laser cooling tech- 
niques that the first Bose-Einstein condensates of atoms 
were achieved 0, 0]. At present, experiments trying to 
achieve ultracold and degenerate molecular gases are cre- 
ating considerable excitement [1, 0, Q. These experi- 
ments should lead to the creation of long-lived molecular 
Feshbach-Einstein condensates, with applications in the 
fields of precision measurements and quantum informa- 
tion 0. 

Near a Feshbach resonance, molecules are formed from 
atoms by tuning a magnetic field and bring into reso- 
nance scattering states of atoms with molecular bound 
states [§]. In this way, conversions between atoms and 
diatomic molecules are induced. A model Hamiltonian 
that describes mixtures of atoms and molecules was in- 
troduced and studied before [^, [l^, Eli [i2| ■ In this let- 
ter, we study this model and analyse the presence or not 
of atomic and/or molecular condensates, using the re- 
cently developed Stochastic Green Function (SGF) algo- 
rithm [13] with tunable directionality IJ]. With this ex- 
act quantum Monte Carlo (QMC) algorithm, momentum 
distribution functions, which are the main indicators of 
condensation, are easily accessible and allow direct com- 
parisons with experiments. We critically compare our 
results with the predictions of mean-field (MF) studies 

13, El. 

We consider the model for bosonic atoms and molecules 
on a lattice. The particles can hop onto neighboring sites, 
and their interactions are described by intra-species and 
inter-species onsite potentials. An additional conversion 
term allows two atoms to turn into a molecule, and vice 



versa. This leads us to consider the Hamiltonian Ti. — 
f + P + C with 

f = -taY^ {a\a^+H.c)-tra ^ {m\m.^+H.c), (1) 
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The T, P, and C operators correspond respectively to 
the kinetic, potential, and conversion energies. The a| 
and Oj operators {to| and m^) are the creation and an- 
nihilation operators of atoms (molecules) on site i, and 
fif = ajoj (n™ = m|mj) counts the number of atoms 
(molecules) on site i. Those operators satisfy the usual 
bosonic commutation rules. The sums {i,j} run over 
pairs of nearest-neighbor sites i and j. We restrict our 
study to one dimension (ID) and we choose the atomic 
hopping parameter = 1 in order to set the energy 
scale, and the molecular hopping parameter t^ = 1/2 
[l^ . The parameter D corresponds to the so-called de- 
tuning in Feshbach resonance physics. In this paper, we 
will systematically use the same value U for the onsite 
repulsion parameters and the conversion parameter. 
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in order to simplify our study. It is important to note 
that the Hamitonian does not conserve the number of 
atoms Na = J2i 'A^^i > nor the number of molecules Nm = 
mlm^, because of the conversion term However 
we can define the total number of particles, N — Na + 
2Nm, which is conserved. 
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While our Hamiltonian is highly non-trivial, it has be- 
come possible to simulate it exactly by making use of the 
SGF algorithm In this algorithm a Green operator 
is considered, 
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where g^q is an optimization matrix, and A are nor- 
malized creation and annihilation operators, defined as 
the operators that create and destroy particles without 
changing the norm of the state they are applied to, and 
{*p|j?} represents two subsets of site indices in which all 
ip are different from all jp. The Green operator is used 
to sample an extended partition function, 

Tr e-^P--'^^ge--^, 
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by propagating across the operator string obtained by 
expanding the exponentials of expression ([6]) in the in- 
teraction picture. When configurations in which Q acts 
as an identity operator occur, then ([6]) reduces to the 
partition function Z{l3) — Tr e^^'^, and measurements 
of physical quantities can be performed. In addition, the 
directionality of the propagation of the Green operator is 
tunable which improves considerably the efficiency 
of the algorithm. 

An important property of the SGF algorithm is that it 
works in the canonical ensemble, the canonical constraint 
being imposed on the total number of particles N. This 
is essential for the efficiency of the simulations. Indeed, 
our model describes a mixture of two different species of 
particles, and would require two different chemical poten- 
tials for a description in the grand-canonical ensemble. 
Adjusting numerically two chemical potentials is cum- 
bersome, because the number of particles of each species 
depends on all parameters of the Hamiltonian. Working 
in the canonical ensemble allows to set the total number 
of particles, and the ratio between the number of atoms 
and molecules is controlled via the detuning, mimicking 
what is done in experiments. 

In order to characterize the different phases encoun- 
tered, it is useful to consider the superfluid density ps- 
The SGF algorithm samples the winding number W, so 
the superfluid density is simply given by ps — {W'^)L/2(3. 
It turns out to be more efficient to measure ps by using 
an improved estimator Wext [13 fo'" the winding number. 
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where are the imaginary time indices of the in- 
teractions appearing when expanding the partition 
function ([5]), and P(Tfe) equals 1 (—1) if a particle jumps 



to the right (left) at imaginary time Tt- This improved 
estimator converges faster to the zero temperature value 
of the winding number [l3|. In our case we evaluate the 
atomic and molecular winding numbers, Wa and Wmi 
and the corresponding superfluid densities are given by 
= (Wl)L/2(3 and = {W^)L/2p. In addition, it 
is useful to define the correlated superfluid density 
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These quantities allow to identify superfluid (SF) and 
super- Mott (SM) [13] phases. SF phases are character- 
ized by non-zero values for p", p™, and pi™ , while a van- 
ishing value for pf'^ with non-zero values for pi and p^ is 
the signature of a SM phase (see caption of FiglT]). The 
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FIG. 1: (Color online) Typical configuration in a SM phase. 
The addition of an extra particle (atom or molecule) has a 
finite energy cost because it creates either a triplet of atoms, 
or an atom-atom-molecule triplet, or a pair of molecules, or 
an atom-molecule pair. Thus the phase is incompressible. 
However exchanging a pair of atoms with a molecule is free, 
thus allowing anti-correlated supercurrents. This correspond 
to third order processes, so the associated superfluid signals 
are small. 



SGF algorithm allows to measure the atomic and molecu- 
lar Green functions (aja^) and (mjmj), from which the 
associated momentum distribution functions na{k) and 
nm{k) are computed by performing a Fourier transfor- 
mation: 

na{k) = — ^ 
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Because we are considering ID systems, we can ex- 
pect at most quasi-condensates. These are characterized 
by a diverging occupation of the zero momentum state 
n{k = 0) as a function of the size L of the system, while 
the condensate fraction n{0)/N vanishes in the thermo- 
dynamic limit. As a result, knowing the value of the con- 
densate fraction for an arbitrary large system size is not 
sufficient to determine if the system is quasi-condensed 
or not. One needs to perform a finite-size scaling anal- 
ysis in order to determine if n(0) diverges or not. In 
the following, all denoted "condensate" phases are to 
be understood as "quasi-condensate" phases. As a re- 
sult, the quantities Uaik) and nm{k) allow us to identify 
phases with atomic condensate (AC), molecular conden- 
sate (MC), or atomic-|-molecular condensates (AC-I-MC). 



3 




FIG. 2: (Color online) Identification of the SM phase. The 
occupation of the zero- momentum states na(0) and nm(0) do 
not change as the size L of the system increases, thus there 
is neither an atomic nor a molecular condensate. However 
the inset shows that the correlated superfluid density pf"^ 
extrapolates to zero as the size increases, while the atomic 
and molecular superfluid densities and con verg e to a 
finite value. This is the signature of a SM phase [13 ]. The 
errorbars are smaller than the L = 60 symbol size. 



It is also useful to keep in mind that the areas below the 
curves na(k) and nm(k) are exactly equal to Na and Nm 
respectively, thus allowing an evaluation of the popula- 
tion of atoms and molecules. 

We concentrate our study on systems with a total den- 
sity ptot — N/L — 2, which is one of the cases con- 
sidered in MF [ni and QMC [H] studies. We investi- 
gate the phase diagram in the {1/U,D) plane. For suffi- 
ciently large interactions U, depending on the detuning 
_D, we find an insulating phase characterized by a vanish- 
ing compressibility and the absence of global superflow. 
This is in agreement with MF studies. However, Ref. 11| 
denotes this insulating phase as a regular Mott insulator 
(MI), whereas we find that it is actually a more exotic 
SM phase (see above). This can be seen in Fig[2]for the 
case U = 1.5 and D ~ —6. The momentum distribution 
functions for atoms and molecules are plotted for differ- 
ent sizes of the lattice. No divergence of the occupation 
of the zero momentum state is perceptible, so there is 
neither an atomic nor a molecular condensate. Moreover, 
the inset shows that the correlated superfluid density p1°^ 
vanishes as the system size increases, as expected for an 
insulating phase. However we can see that the super- 
fluid densities associated with the individual atomic and 
molecular species converge to a finite value, which is the 
signature of a SM phase (a similar phase is also present 
in the case of Bose-Fermi mixtures [l5| ) . This is the first 
qualitative difference between MF results and ours. 

Starting from the above SM phase, reducing the in- 
teractions will eventually break the solid structure. For 
U = \ and negative detuning _D = — 6 the system under- 
goes a molecular condensation, as can be seen in FiglH] 
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FIG. 3: (Color online) Identification of the MC phase. The 
occupation of the zero- momentum state of atoms UaiQ) does 
not change as the size L of the system increases, thus there 
is no atomic condensate. However the molecular momentum 
distribution function becomes narrow and nm(0) diverges as 
the size of the system increases, which is the signature of a 
MC phase. The inset shows the L = 60 results separately 
for clarity. The errorbars are smaller than the L = 60 symbol 
size, except in fe = for molecules for which the error is about 
2 times the size of the symbol. 



No divergence of the occupation of the zero momentum 
state na(0) occurs. However nm(0) diverges and shows 
the presence of a molecular condensate. This transition 
from an insulator to a MC phase is in agreement with 
MF theory. 

The SM phase persists when going from negative to 
positive detuning with large interactions. For sufSciently 
large detuning D, MF studies predict a direct transition 
from a MI (actually SM) phase to an AC-I-MC phase, 
as the interactions are reduced. However, our remark- 
able result is that we find an intermediate AC phase in 
a small region of the phase diagram. This can be seen 
in Fig [4] for the case [/ = 5 and D — A. We can see 
that a divergence of ^^(O) occurs while nm(0) remains 
constant as the system size increases. Thus we are in 
presence of a phase in which the atoms are condensed, 
but not the molecules. Such an AC phase is missing, to 
our knowledge, in all phase diagrams coming from MF 
theory 1^, 11 1 . Our evidence for the AC phase in Fig. [4] 
is comparable to that for the MC phase in Fig. [31 

While it is hard to give a phase diagram with pre- 
cise borders delimiting the different phases (because each 
point of the diagram requires a heavy finite-size scaling 
analysis), we give a qualitative diagram in FigObased on 
simulations for considerably more values of the detuning 
D and the interaction U than in Figl2| to [H In addition, 
we provide connection with future experiments by show- 
ing on Figinihow the atomic and molecular visibilities Va 
and Vm ^16] behave when entering the AC phase from the 
SM phase along the vertical line D = —2. In the present 
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FIG. 4: (Color online) Identification of the AC phase. The 
occupation of the zero-momentum state of molecules nm{0) 
does not change as the size L of the system increases, thus 
there is no molecular condensate. However the atomic mo- 
mentum distribution function becomes narrow and na{0) di- 
verges as the size of the system increases, which is the signa- 
ture of an AC phase. The inset shows the L = 60 results sep- 
arately for clarity. The errorbars are smaller than the L = 60 
symbol size, except in = for atoms for which the error is 
about 4 times the size of the symbol. 



case Va reaches unity before Vm as the interactions de- 
crease, showing again that the system is entering an AC 
phase. 

To conclude, we have studied a two-species Bose- 
Hubbard Hamiltonian for atoms and molecules on a lat- 
tice, interacting via a Feshbach resonance. We have 
shown that the MI phase identified in MF studies is ac- 
tually a SM phase. For large negative detuning we find 
a transition from this insulating phase to a MC phase, 
in agreement with MF theory. For smaller negative or 




FIG. 6: (Color online) The atomic and molecular visibilities 
Va and Vm- Because of the divergence of na{0) or nm{0), Va 
(and/or Vm) must converge to unity when an atomic (and/or 
a molecular) condensate occurs. 



positive detuning however, while MF theory predicts a 
direct transition from MI to AC-I-MC, we find that an 
AC phase occurs and that the system undergoes phase 
transitions from SM to AC to AC+MC. The AC phase 
was not found in previous studies. The phase diagram 
we provide may serve as a guide for the investigation of 
atomic and molecular quantum matter now that conden- 
sation of Feshbach molecules is begining to be achieved 
experimentally. 
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FIG. 5: (Color online) The qualitative phase diagram in the 
(1/(7, D) plane for ptot = 2. We identify regions with super- 
Mott (SM), atomic condensate (AC), molecular condensate 
(MC), and atomic -l-molecular condensate (AC-I-MC) phases. 
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